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1 Introduction 



The internal structure of the nucleon is the most fundamental problem of strong 
interaction physics. Understanding this structure in terms of the elementary 
quark and gluon degrees of freedom of the underlying theory, quantum chro- 
modynamics (QCD), remains the greatest unsolved problem of the Standard 
Model of nuclear and particle physics. 

Historically, the basic strong interaction which we have sought to explain 
has been that between protons and neutrons in the atomic nucleus. The origi- 
nal idea of massive particle exchange of Yukawata has been a guiding principle 
according to which later theories have been developed. It was pointed out by 
Wick □ that this idea was consistent with the Heisenberg Uncertainty princi- 
ple, whereby the interaction range of the nuclear force is inversely proportional 
to the mass of the exchanged meson. Over the years a phenomenological de- 
scription of the forces acting between nucleons has been developed within a 
meson-exchange picture. 

Following the experimental discovery of the pion in 1947, the 1950s and 
1960s saw an explosion of newly discovered mesons and baryons, as particle 
accelerators pushed to higher energies. To bring some sense of order to the pro- 
fusion of new particles, Gell-Mann and Zweig introduced the idea of quarkstJ, 
which enabled much of the hadronic spectrum to be organized in terms of just 
a few elementary constituents. Soon after, however, it was realized that a se- 
rious problem existed with the simple quark classifications, namely the A ++ 
isobar. The quark model wave function for the A ++ was predicted to be totally 
symmetric, however the A ++ obeyed Fermi-Dirac statistics. A solution to this 
problem was found by assigning extra internal colora quantum numbers to the 
quarks, in which baryons would have in addition an antisymmetric color wave 
function. 

The discovery of scaling in deep-inelastic electron-nucleon scattering in 
the late 1960s at SLACH confirmed that the nucleon contained point-like con- 
stituents, which were soon identified with the quarks of the quark model. Im- 
posing local gauge invariance on the color fields, and introducing vector gluon 
exchange to mediate the inter-quark interaction, led naturally to the develop- 
ment of QCD as the fundamental theory of strong interactions 

Because QCD is an asymptotically free theory — the effective strong cou- 
pling constant decreases at short distances — processes involving large mo- 
mentum transfers can be calculated reliably within perturbation theory. Yet 
despite the successes of perturbative QCD, we are still unable to extract from 
QCD sufficient details regarding its long-distance properties. This is because 
in the infra-red region the strong coupling constant grows and perturbation 
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theory breaks down, and the available non-perturbative tools are not yet suf- 
ficiently developed to allow quantitative predictions. 

In a sense it is ironic that the theory which arose out of the desire to 
understand nuclear forces is able to explain backgrounds in hadronic jets pro- 
duced in high energy collisions, yet is unable to describe the properties of the 
ground state of the theory. Although one can argue that QCD in principle 
explains all hadronic and nuclear phenomena, without understanding the con- 
sequences of QCD for hadron phenomenology one may as well argue that the 
entire physics of atoms and molecules can in principle be explained from QEDI3. 
Understanding how the transition from the quarks and gluons of QCD to the 
physical mesons and baryons takes place remains the holy grail of modern 
nuclear physics. 

In the next Section some basic elements of QCD relevant for later appli- 
cation to nuclcon structure are reviewed. This is followed in Section 3 by the 
basic definitions and kinematics of elcctron-nuclcon scattering, including elas- 
tic, deep-inelastic and semi-inclusive scattering. In Section 4 we focus more 
closely on the flavor and spin content of the nucleon, and outline some recent 
highlights in the study of valence and sea quark distributions. Finally, some 
concluding remarks are made in Section 5. 



2 Elements of QCD 

Quantum Chromodynamics is a non-Abelian gauge field theory based on the 
gauge group SU(3) co ior, and defined in terms of the Lagrange density 6 > 8 ~ 10 : 

£qcd = Anv + £ gauge ~r ^ghost i (i) 

where £; nv is the classical Lagrangian, invariant under local gauge transforma- 
tions of the SU(3)coior group: 

£ inv = \ f {i^D^ - m f ) l3 ^j - \f; v F^ a . (2) 

The quark fields iftij (for a particular quark flavor / = u,d, s, ■ ■ ■ , with mass 
ruf) are labeled by color indices i,j = 1,2,3. The covariant derivative is 
Dfi = 9ju — ig^A't, where g is the QCD coupling constant and T a are the 
generators of the SU(3) group, with a = 1, • • • , 8. In terms of the gluon field 
the gluon field strength tensor is F£ v = d^A a v - d v A% + gj ahc A h ^A c vl with 
f abc the SU(3) structure constants. The major differences between QCD and 
quantum electrodynamics is the appearance, due to the non-Abelian structure 
of the theory, of gluon self-couplings in the F ■ F term. This gives rise to 3- 
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and 4-point gluon interactions which make the theory highly non-linear, but 
also leads to the property of asymptotic freedom (Section 3.2). 

Under local gauge rotations the quark fields transform (dropping color and 
flavor indices) according to: 

r/>(x) -> tf(x) = U(x) 1>{x) , (3) 

where U is an SU(3) unitary matrix: 

U(x) = exp{id a (x)T a ) , (4) 

with a (x) real. The gluon fields transform according to: 

A M (x) -» A'^x) = U(x)A ll (x)U- 1 (x) + I (d^U(x)) CT^x) . (5) 

One can easily show that C- mv is then invariant under these transformations. 

The gauge invariance of the classical Lagrangian introduces some diffi- 
culties when quantizing the gauge theory. This problem is avoided with the 
introduction of an addition term, £ ga ugc, which fixes a specific gauge. In 
the Lorentz (covariant) gauge, one has £ ga ugo = —(l/2a)(d l *A*) 2 , where a is 
an arbitrary gauge parameter. Of course observables cannot depend on the 
choice of a, and some common choices are a = 1 (Feynman gauge) and a — > 
(Landau gauge). Other, non-covariant, gauge choices are the Coulomb gauge 
{diAf = 0), the axial gauge (A3 = 0), and the temporal gauge (Ag = 0). 

The Faddeev-Popov ghost density, £ ghost = (<9 M x a )(c5 afc <9 M - gf abc A L ^)x c , 
where x a arL d x a are scalar aniz-commuting ghost fields, ensures that the gauge 
fixing does not spoil the unitarity of the S-matux- Further discussion about 
the gauge fixing problem can be found in Refs.liiM. 

In renormalizable field theories such as QCD the strength of the interaction 
depends on the energy scale. A property almost unique to QCD is that the 
renormalized coupling constant decreases with energy — known as asymptotic 
freedom. The running of the QCD coupling with energy allows one to compute 
cross sections for any quark-gluon process using a perturbative expansion if 
the coupling is small. Writing the full QCD Lagrangian as £qcd = £0 + A) 
where 
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Co = ^{i^d^ - m)V - ^(MS - d u A a ^)(d' i A a v - d v A a ") 

-^(^) 2 + (c^x a )(^x a ) (6) 
is the free Lagrangian, and 

— ^f abe f cde A a ^A b J A cll A b ' J - gf abc (d^x a )x b ^ (7) 

the interaction part, one can derive from C\ a complete set of Feynman rules 
for computing any scattering amplitude involving quarks and gluons 8-11 . 

Perturbative QCD has been enormously successful in calculating hard pro- 
cesses in high energy lepton-lepton, lepton-hadron and hadron-hadron scatter- 
ing. However, even at high energies one can never avoid the fact that the 
physical states from which the quarks and gluons emerge to undergo the hard 
scattering are hadrons, so that one always encounters soft scales in any strongly 
interacting system. While operator product expansions usually allow one to 
factorize the short and long distance dynamics, understanding the complete 
physical process necessarily requires going beyond perturbation theory. 

Over the years the problem of non-perturbative QCD has been tackled on 
several fronts. The most direct way is to solvepthe QCD equations of motion 
numerically on a discretized space-time lattice 113. Recent advances in lattice 
gauge field theory and computing power has made quantitative comparison of 
full lattice QCD calculations with observables within reach. 

Alternative methods of tackling non-perturbative QCD involve the build- 
ing of soluble, low-energy QCD-inspired models, which incorporate some, but 
not all, of the elements of QCD. Phenomenological input is then used to con- 
strain the model parameters, and identify circumstances where various approx- 
imations may be appropriate. These approaches often exploit specific symme- 
tries of QCD, which for some observables may bring out the essential aspects 
of the physics independent of the approximations used elsewhere. A good ex- 
ample of this, which has had extensive applications in low energy physics, is 
chiral symmetry. 

Consider the classical quark Lagrange density in Eq. (|^) in the limit where 
the mass of the quarks is zero: 
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£L = W-fD^ = i> L i^D^ L + ^Hi-fD^R , (8) 

where ipL.R = (1/2)(1 ± 75)1/' are left- and right-handed projections of the 
Dirac fields. Under independent global left- and right-handed rotations Cf nv 
remains unchanged. For Nf massless quarks, the classical QCD Lagrangian is 
then said to have a chiral SU(iVf )i<g> SU (Nf)n symmetry. Of course non-zero 
quark masses break this symmetry explicitly by mixing left- and right-handed 
quark fields, however, for the u and d quarks, and to some extent the s, the 
masses are small enough for the chiral symmetry to be approximately valid. 

If chiral symmetry were exact, a natural consequence would be parity 
doubling. The nucleon would have a negative parity partner with the same 
mass, and the pseudoscalar pion would have the same mass as the scalar meson. 
In nature, the lightest negative parity spin- 1/2 baryon is the Six resonance, 
which is several hundred MeV heavier than the nucleon. Moreover, the pion has 
an exceptionally small mass, while the lightest candidate for a scalar meson 
is several times heavier than the pion, so that chiral symmetry in nature is 
clearly broken. 

The way to reconcile a symmetry which is respected by the Lagrangian but 
broken by the physical ground state is if the symmetry is broken spontaneously. 
According to Goldstone's theorem, a consequence of a spontaneously broken 
chiral symmetry is the appearance of massless pseudoscalar bosons. For Nf = 2 
(namely, for u and d flavors), these correspond to the pseudoscalar pions; for 
Nf — 3 (counting the strange quark as light), these include in addition the 
kaons and the 77 meson. The physical mesons are of course not massless, but 
on the scale of typical hadronic masses (~ 1 GeV), they can be considered 
light — their masses arising from the small but non-zero quark masses. A 
perturbative expansion in terms of the small pseudoscalar boson masses can 
be developed Ej, and applied systematically to describe hadron interactions at 
very low energies. 

As will be demonstrated in Sections 4.2 and 4.3 the chiral properties of 
QCD are in fact critical to understanding many aspects of nucleon structure, 
from low energy form factors to deep-inelastic structure functions. Before pro- 
ceedings with further discussion about QCD and nucleon structure, however, 
in the next Section we first define the observables through which one can study 
the internal structure of the nucleon in electron scattering. 
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3 Electron Nucleon Scattering 

Because the electromagnetic interaction of leptons is perhaps the best un- 
derstood part of the Standard Model, the cleanest way to probe the internal 
structure of hadrons is through lepton scattering. This applies to both charged 
leptons and neutrinos, although to be concrete we shall consider the scattering 
of electrons. 

In the one-photon approximation, the scattering of an electron with four- 
momentum I from a nucleon with momentum p is illustrated in Fig. 1, where 
the outgoing electron momentum is V and the hadronic final state is denoted 
by X: eN — > e'X. The energies of the incident and scattered electrons are 
E and E', and the electron scattering angle is 9. The energy transfer to the 
nucleon in the target rest frame is v = E—E', and the four- momentum transfer 
is q 2 = (l- I') 2 ~ -4EE' sin 2 6 for m e <^E,E'. 

In the nucleon rest frame the differential cross section is given by: 



d 2 a a 2 ern E' 



dttdE' 2MQ 4 E 



TV , ( 9 ) 



where a em is the electromagnetic fine structure constant, M is the nucleon 
mass, and Q 2 = — q 2 . The normalization of states in Eq.(^) is such that 
(p\P ') = (2tt) 3 2po 5(p — p'). The lepton tensor L^ v is given by: 

jju> = 2 lH iv + 2 + g ^ q 2 T 2ie" uXp l x l' p , (10) 

corresponding to an electron with helicity ±1/2. Note that for unpolarized 
scattering only the part of L^ v which is symmetric under the interchange 
fi <-> v is relevant, while for polarized only the antisymmetric part enters. 
The hadronic tensor, 

iV = i]C(2^ 4 (p + g -p x ) WM\x)WM\$) > (n) 

x 

contains all of the information on the structure of the target. For inclusive 
scattering one sums over all final states X, while for exclusive scattering X 
denotes a specific hadron. 

The most general form for the hadronic tensor consistent with Lorentz and 
gauge invariance, as well as invariance under time reversal and parity, is: 

^W, v = {-g, v + Q -f-)w 1+ (p, ^) (p„ - P -^q v ) W 2 

+ ie^X P jj(s p M 2 G 1 + (p-q S p -s-qpP)G 2 ) , (12) 
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Figure 1: Electron-nucleon scattering in the one-photon exchange approximation: X = N 
for elastic scattering, while X is a sum over hadrons for inclusive inelastic scattering. 



where s p is the nucleon spin vector, defined by s p (X) — (2X/M)(\p \;poP) for 
nucleon helicity A, so that s 2 — — 1, s -p = 0. The structure functions W%, W 2 , 
Gi and G 2 are in general functions of two variables. Usually one chooses these 
to be Q 2 , and the Bjorken x variable, defined as x = Q 2 /2p ■ q. 

Since W± and W% are coefficients of Lorentz tensors symmetric in /iv, they 
can be measured in unpolarized electron-nucleon scattering. The unpolarized 
differential cross section can be written in the target rest frame astfl'lij: 

j2 Ttr+ltr sn,2 pi? 
dndE , = -qT- [Wx^Q^sin^m + W^Q 2 )] , (13) 

where ttl" (lit) refers to the polarization of the electron parallel (antiparallel) 
to that of the target nucleon. 

The structure functions G\ and Gi can be measured by taking the dif- 
ference of ccass. sections with electron and nucleon polarizations parallel and 



[{E + E'cos6)MG 1 {x,Q 2 )-Q 2 G 2 (x,Q 2 )] . (14) 



4 15 



antiparallel I 

dVLdE 1 Q 4 E 



It will be convenient later to introduce dimensionless structure functions 
F\,2, 5i,2, defined as: 

M Wi = Fi, vW 2 = F 2 

M 2 v Gi = si, Mv 2 G 2 = 92 ■ (15) 

As we will see in Section 3.2, some of these dimensionless structure functions 
have very simple interpretations in terms of quark densities in deep-inelastic 
scattering. 
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Figure 2: Electron— nucleon differential cross section as a function of hadronic final state 



The structure functions describe scattering to final states whose spectrum 
depends on the amount of energy and momentum transferred to the target 
nucleon. In Fig. 2 the differential cross section for unpolarized scattering is 
plotted as a function of the invariant mass, W, of the hadronic final state, where 
W 2 = (p + q) 2 = M 2 + 2Mu - Q 2 . At low energy, since v > Q 2 /2M only 
elastic scattering is kinematically allowed (represented by the spike in Fig. 2 at 
W = M). As the energy increases above the pion production threshold, Wth = 
M +m 7 r (left vertical line in Fig. 2), inelastic scattering to nucleon + multi-pion 
states can occur, as well as excitation of nucleon resonances. The first peak 
corresponds to the spin- and isospin-3/2 A resonance at W = 1232 GeV. The 
second peak is predominantly due to the negative parity partner of the nucleon, 
the Sn resonance, and the third to the Fi 5 . In the region between W = W t h 
and W — 2 GeV many resonances contribute, most of whose contributions are 
buried underneath the background. The vertical line in Fig. 2 at W — 2 GeV 
corresponds to the approximate boundary between the resonance and deep- 
inelastic scattering (DIS) regions. In the next Section we shall examine electron 
scattering to the simplest final state, namely the elastic. 

3.1 Elastic Form Factors 

The most basic observables which reflect the composite nature of the nucleon 
are its electromagnetic form factors. Historically, the first indication that the 
nucleon is not elementary came _£com measurements of the form factors in 



mass, W . 



elastic electron-proton scattering! 
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The nucleon form factors are defined through matrix elements of the elec- 
tromagnetic current, J M = ipj^tp, with the quark field, as: 

{N{P')\J^)\N{P)) = u(P') ( 7M *i (Q 2 ) + l -^fT 2 {Q 2 )^ U {P) , (16) 

where P and P' are the initial and final nucleon momenta, and q = P — P' . 
The functions T\ and T 2 are the Dirac and Pauli form factors, respectively. In 
terms of T\ and Ti the Sachs electric and magnetic form factors are defined 
as: 

G E (Q 2 ) = ^i(Q 2 ) - (Q 2 /4M 2 ) T 2 {Q 2 ) , (17) 
G M (Q 2 ) =^i(Q 2 )+^ 2 (Q 2 ) • (18) 

Squaring the amplitude in Eq.(|l^) and comparing with the cross section in 
Eq. ( |l3| ) , one can write the structure functions for elastic scattering in terms of 
the electromagnetic form factors: 

i2 fr,2 



F£ = Mt G M (Q ) 8 \ v — J , (19) 

= (G 2 e {Q 2 ) + rGlAQ 2 )) 6 (v - , (20) 

where r = Q 2 /AM 2 . For the spin-dependent structure functions one has: 

fjf = ^G M (Q 2 )(G E (Q 2 )+rG M (Q 2 ))S^-^j , (21) 

df = ^G M (Q 2 ) (G E {Q 2 ) - G M (Q 2 )) s(u-^j . (22) 

At the quark level, the form factors can be decomposed as: 

G e , m (Q 2 )^Y,^ G e.,m(Q 2 ) > ( 23 ) 

so that contributions from specific quark flavors can be identified by considering 
different hadrons. For the proton one has: 

2 1 1 

f<P _ fin fid f-is (OA\ 

while for the neutron (using isospin symmetry): 

fin ^ fid 1 fiu ^ fis {OK\ 
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where G q E M by definition refers to the quark flavor in the proton. 

At low Q 2 the distance scales on which the electromagnetic scattering takes 
place are large, so that the resolving power of the probe is only sufficient to 
measure the static properties of the nucleon. In this region the form factors 
reflect completely the non-perturbative, long-distance structure of the nucleon. 
At Q 2 = the electric form factor is equal to the charge of the nucleon, while 
the magnetic form factor gives the magnetic moment: 

G£(0) - e N , G£(0) - m , (26) 

where = 1(0) for the proton (neutron) , and in units of the Bohr magneton, 
p p = 2.79^o an d M« = — 1-91/io- 

At small Q 2 away from zero, in a frame of reference where the energy 
transfer to the nucleon is zero (namely, the Breit frame, v = 0, Q 2 = q 2 ), 
the electric and magnetic form factors measure the Fourier transforms of the 
distributions of charge and magnetization, pe and pm, in the nucleon: 

G e ,m(Q 2 ) = J d 3 r e-^ p PE . M (r) . (27) 

Expanding the form factors about Q 2 — 0, 

G e ,m(Q 2 ) = G e ,m(0) - \q 2 {t 2 )e,m + G{Q A ) , (28) 

6 

enables one to define the charge and magnetization radii of the nucleon in terms 
of the slope of the form factors at Q 2 = 0: (t 2 )e,m = —6 (1Ge,m /dQ 2 \ Q2 _ Q . 
Empirically, the form factors at low Q 2 are well described by a dipole form, 




where Qq = 0.71 GcV 2 , in which case the electric and magnetic r.m.s. radii 
of the proton, and the magnetic radius of the neutron, are (r 2 ) 1 / 2 = 0.81 fm. 
Because the neutron has zero charge, the neutron electric form factor, although 
non-zero, is very small. 

The dipole form can be qualitatively understood within a vector meson 
dominance picture, in which the photon at low Q 2 fluctuates into a J p = 
1~ meson (such as the p), which then interacts with the nucleon. The Q 2 
dependence of the form factor is then given by the vector meson propagator, 
with the mass of the p corresponding approximately to Qq. However, while 
providing a reasonable approximation to the form factors at low Q 2 , deviations 
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from the dipole form have been observed, and it is important to understand 
the nature of the deviations at larger Q 2 . 

At the other extreme of asymptotically largaQ 2 , the elastic form factors 
can be described in terms of perturbative QCDE3. Here the short wavelength 
of the highly virtual photon enables the quark substructure of the nucleon 
to be cleanly resolved. By counting the minimal number of hard gluons ex- 
changed between quarks, the Q 2 behavior of the form factors is predicted from 
perturbative QCD to be Ge,m{Q 2 ) ~ 1/Q 4 at large Q 2 . Just where the per- 
turbative behavior sets in, however, is still an open question which must be 
resolved experimentally Evidence from recent experiments at Jefferson Lab 
and elsewhere suggests that non-perturbative effects are still quite important 
for Q 2 at least » 10 GeV 2 . 

Understanding the transition from the low to high Q 2 regions is vital 
not only for determining the onset of perturbative behavior. Form factors in 
the transition region at intermediate Q 2 are very sensitive to mechanisms of 
spin-flavor symmetry breaking, which cannot be described within perturba- 
tion theory. In Section 4 we give several examples where form factors reflect 
important aspects of the non-perturbative structure of the nucleon. 

3.2 Deep -Inelastic Structure Functions 

Because of the 1/Q 4 dependence in the elastic form factor in Eq. (p9|) , the 
elastic cross section dies out very rapidly with Q 2 . It was therefore expected 
that the inelastic cross section would behave in a similar fashion at large Q 2 . 
Contrary to the expectation, the observational in the late 1960s that the in- 
elastic structure function does not vanish with Q 2 , but remains approximately 
constant beyond Q 2 ~ 2_GeV 2 , provided the first evidence of point- like con- 
stituents of the nucleon PHr ! and led to the development of the parton model 
and later QCD. 

For inclusive scattering the sum in Eq. ( |ll| ) is taken over all hadronic final 
states X. Using the completeness relation Ylx I J ^)( J ^I = 1 an d translational 
invariance, the hadronic tensor can be written: 

w » v = hf d ^ eiq< (mMOM0)\N) , (30) 

where £ is the space-time coordinate. In the limit where p-q and Q 2 — > oo, but 
the ratio of these fixed (Bjorkcn limit), receives its dominant contributions 
from the light-cone region. This is clear if one writes the argument of the 
exponential in light-cone coordinates, 

9 • f = ~Y~ + ~~ 2 1-L " d- i (31) 
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where £± = ^ ± In the target rest frame the photon momentum can 
be taken as = (v; 0±, — \Jv 2 + Q 2 ^j ~ (v; 0^, — v — Mx), so that g • £ = 

-Ma;(£o - 6*)/2 + {2v + Mx)(£ + £ z )/2, where in this frame x = Q 2 /2Mv. 
The largest contributions to the integral are those for which the exponent 
oscillates least, namely q ■ £ ~ 0. In the Bjorken limit q ■ £ behaves like i/(£o + 
£ z ), so that only when £o = will there be non- negligible contributions to 
W^v Therefore the DIS cross sections is controlled by the product of currents 
Jn(0 Jv(0) near the light cone, £ 2 ~ 0. 

Parton Model 

The connection between the deep-inelastic structure functions and the_quark 
structure of the nucleon was first provided by Feynman's parton modelcJ. The 
hypothesis of the parton model is that the inelastic scattering is described by 
incoherent elastic scattering from point-like, spin 1/2 constituents (partons) in 
the nucleon. 

The validity of the parton picture relies on the treatment of the interactions 
of the virtual photon with the partons in the impulse approximation. The 
legitimacy of the impulse approximation rests on two assumptions: (i) final 
state interactions can be neglected, and (ii) the interaction time is less than the 
lifetime of the virtual state of the nucleon as a sum of its on-shell constituents. 

The first assumption seems reasonable since in DIS the energy trans- 
ferred to the parton is much greater than the binding energy, so that the 
partons can be viewed as quasi-free. The second assumption can be verified 
in the infinite momentum frame (IMF), where the momentum of the nucleon 
is = (p z + M 2 /2p z ; 0j_, —p z ), with p z — > oo (or v/c — > 1), in which case 
the photon four-momentum is q^ = (— xp z (l — M 2 x 2 /Q 2 ) + Mv/2p z \ 0^, 
xp z (l — M 2 x 2 /Q 2 ) + Mv/2p z J. Since the dominant contributions to W^v are 
those for which q ■ £ -c 1 and £o — — £z, one has q ■ £ ~ 2Mv£o/p z , so that the 
interaction time is ~ £o < p z /2Mv. The lifetime of the virtual state can be 
obtained by observing that the energy of the virtual nucleon state consisting of 
on-shell partons with momenta Xip z and mass mi is ~ ^2ii x iPz + Trq /2xip z ), 
so that the difference between the energies of the virtual and on-shell nuclc- 
ons is s» (J2i m i l x i ~ M 2 )/2Pz- Therefore the lifetime of this virtual state is 
proportional to p z , and in the Bjorken limit the ratio of interaction time to 
virtual state lifetime ~ 1/v — > 0. 

The parton picture is then of a quark with momentum fraction n absorbing 
a photon with Xi = x, since S((q + Xip) 2 ) — > 8(x — Xi)/2p ■ q. One can then 
relate the structure functions to the parton densities (quark and antiquark 
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momentum distribution functions) in the nucleon as 



0@ 



F 2 (x) = 4 x(q(x)+q(x)) = 2xF 1 (x) , 



(32) 



9l{x) = ^Y, e l x ( A <i( x ) + A ^)) > 



(33) 



where q(x) = q*(x) + q-^(x) and Aq(x) = q*(x) — q^(x) are the spin-averaged 
and spin-dependent quark densities. The consequence of point-like partons 
is the non-vanishing of the inelastic structure functions at large momentum 
transfers, since the structure functions in Eqs.(|32|) and (|3^) are independent 



Although providing a simple, intuitive language in which to interpret the 
qualitative features of the deep-inelastic data, the parton model is not a field 
theory. The formal basis for the parton model is provided by the operator 
product expansion and the renormalization group equations in QCD, which 
actually gives rise to small violations of scaling through QCD radiative correc- 
tions. 

Operator Product Expansion 

In quantum field theory products of operators at the same space-time point 
(composite operators) are not well defined lid. The short distance operator 
product expansion (OPE) of Wilson Ell, in which the composite operators are 
expanded in a series of finite local operators multiplied by singular coefficient 
functions, provides a way of obtaining meaningful results. 

Because deep-inelastic scattering probes the £ 2 ~ region, rather than 
the £ ~ 0, one needs an expansion of the product of currents in Eq.(|30|) which 
is valid near the light-cone (this is because at short distances Q — > oo and 
p ■ q/Q 2 — > 0, while in DIS the light-cone region corresponds to the Bjorken 
limit, Q 2 — » oo and p ■ q/Q 2 ^ 0(1)). The general form of the light-cone 
operator product expansion isEH: 



where the sum is over different types of operators with spin N (i.e. those what 
transform as tensors of rank iV under Lorentz transformations). In DIS the 
spin-iV operators O^ 1 MJV represent the soft, non-perturbative, physics, while 
the coefficient functions describe the hard photon-quark interaction, and 
are calculable within perturbative QCD. 



of Q 2 . 



(34) 



i,N 
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It is useful to categorize the operators according to their flavor properties, 
namely those that are invariant under SU(iV/) flavor transformations (singlet) 
and those that are not (non-singlet). For unpolarized scattering (the exten- 
sion to spin-dependent scattering is straightforward) the operators must be 
completely symmetric with respect to interchange of indices^ i ■ • • /xjvj so that 
one can construct at most 3 kinds of composite operatorsEH13. The non-singlet 
operators must be bilinear in the quark fields: 

O^""" = j-j^j ip ( 7 ^ 1 D M2 • • • D» N + fiifij permutations) A ijj , (35) 

where A are the eight Gell-Mann matrices of the flavor SU(iV/) group. The 
singlet operators are: 

jN-l_ 

_ __^( 7 Mi jD M2 ... d »n + ^ permutations) ip , (36) 

■JV-2 

g'"^ N = Yn\ ( F>liaD ^ ■ • • D^N-ipgN + ^ permutations) , (37) 

corresponding to the quark and gluon fields, respectively, and color indices 
have been suppressed. 

Equations (j35"l)-(|37|) represent operators with the lowest 'twist', defined as 
the difference between the mass dimension and the spin, N, of the operator. 
Whereas the leading twist terms involve freexuiark fields, operators with higher 
twist involve both quark-gluon interactionst3, for example <t/jF liV/ y t/ ij), and are 
suppressed by powers of 1/Q 2 . 

The matrix elements of the operators contain information about the long- 
distance, non-perturbative structure of the nucleon. They can in general be 
written as: 

{N(p)\0? 1 '" ,tN \N(p)) = Af p^ 1 ■ ■ -p» N - (g^ terms) , (38) 

where Af represents the soft physics, and the 'trace terms' containing the 
gf-iH a re necessary to ensure that the matrix elements are traceless (i.e. so 
that the composite operator has definite spin, N). When contracted with the 
these give rise to terms that contain smaller powers of v 2 (i.e. Q 2 = Oiy) 
instead of (p ■ q) 2 = 0{v 2 )). 

The OPE analysis allows one to factorize the moments of the structure 
functions into short and long distance contributions, where the latter are 
target-dependent (and Q 2 -independent). For the F 2 structure function, for 
example, one has: 

M 2 N (Q 2 ) ee fdx x N ~ 2 F 2 (x,Q 2 ) = ^Cf(Q 2 ) Af . (39) 
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The target-independent (and Q 2 -dependent) coefficient functions can be 
calculated at a given order in perturbation theory directly from the renor- 
malization group equations, which will introduce logarithmic Q 2 violations of 
scaling compared with the simple parton model. 

The structure function can be obtained from the moments via the inverse 
Mellin transform: 

i rNn+ioc 

F 2 (x, Q 2 ) = — / dN x^ N M 2 N (Q 2 ) , (40) 

11:1 JN -ioc 

by fixing the contour of integration to lie to the right of all singularities of 
M2 (Q 2 ) in the complex- N plane. In terms of the quark and gluon distribu- 
tions, F2 ja given by a convolution of the parton densities with the coefficient 
functions describing the hard photon-parton interaction: 

F p 2 {x,Q 2 ) = \£ J (c NS (x,Q 2 ) xq NS (^,Q 2 ) 

+ ^ Cq (x,Q 2 ) xS^,Q 2 )+ \c G ( Xl Q 2 ) *g(~,Q 2 )) , (41) 

where qNs(x, Q 2 ) — (u + u — d — d — s — s)(x, Q 2 ) is the flavor non-singlet 
combination (for three flavors), while the singlet combination T,(x,Q 2 ) = 
J2 q (l + q)( x ,Q 2 ): an d G(x,Q 2 ) is the gluon distribution. The gluon coeffi- 
cient Cq enters only at order a s (since the photon can only couple to the 
gluon via a quark loop). 

The challenge to understanding the quark and gluon structure of hadrons 
is to calculate the soft matrix elements .4^ in Eq.(|39|), or equivalently the 
parton distributions in Eq.(f4l|), from QCD. At present this can be only be 
done numerically through lattice QCD, or in QCD-inspired quark models of 
the nucleon. Considerable progress has been made over the last two decades 
in attempting to establish a connection between the high energy parton pic- 
ture of DIS on the one hand, and the valency nuark models at low energy on 
the other 24 ~ 26 . The underlying philosophy EZro has been that if the nucleon 
behaves like three valence quarks at some low momentum scale ~ fi 2 , a purely 
valence quark model may yield reliable twist-two structure functions. Com- 
parison with experiment at DIS scales, where a description in terms of valence 
quarks will no longer be accurate, can then be made by evolving the structure 
function to higher Q 2 via the renormalization group equations. 
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Renormalization Group Equations 

In an interacting field theory like QCD quantities such as coupling constants, 
masses, as well as wave functions (operators), must be renormalized. The 
renormalization procedure introduces an arbitrary renormalization scale, fj?, 
into the theory, although of course the physics itself cannot depend on /i 2 . 

In the following we will consider the renormalization of the non-singlet 
operators corresponding to the F2 structure function. The generalization to 
singlet operators in straightforward, although one needs to take into account 
mixing between the singlet quark and gluon operators£3. If the unrenormalizcd 
matrix elements of the OPE are independent of the renormalization scale [A, 
then 

±(N(p)\Mt)U0)\N(p))=0- (42) 

Defining the wave function renormalization of the spin- TV non-singlet operator 
by 0j^ rc = Z N 0^ n , where Z N is the renormalization constant, Eq. (f42|) can 
be rewritten as: 

(m^ + m-jjj - 7*) C N (Q 2 /n, g) = 0, (43) 

for each spin N. This is the well-known renormalization group equation for the 
coefficient functions. In Eq.(^), the strong coupling constant g is renormalized 
at the scale /i 2 , and ^ N is the anomalous dimension of the twist- two operator 
O n : 

1 "=^±(lnZ N ), (44) 
and the /3-function is given by: 



The solution to Eq.© is 

where g is the effective (running) coupling constant, defined by dg 2 /dt — 
g (3(g), with t = ln(Q 2 //i 2 ) and g(t — 0) = g. The coefficients C N , anoma- 
lous dimension j N (g) and the (3- function can all be calculated in perturbation 



(45) 
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theory by expanding in powers of the coupling, g: 



7 Ar (ff)=7 (0)A '/ I + 0(g 2 ), 

0(9) = -Po + 0(g 5 ) , (47) 

C W (1 ! 5 2 ) = C(°) JV + 0( 3 2 ). 

It is then straightforward to derive the equation governing the Q 2 evolution of 
the .pan-singlet moments of the structure functions, which to lowest order in g 



Q N (Q 2 ) = I' dxx N - 1 q NS (x) 
Jo 



Q N (V 2 ) , (48) 



a s (^ 2 ) 

where the lowest order non-singlet anomalous dimension is: 




(49) 



and /3q = 11 — 2Nf/3 for Nf active flavors in the evolution. The strong coupling 
constant has been rewritten in Eq.(B§) as: 



(n 2, _ f(Q 2 ) 47T 

a s {Q) = — = — , n2/ . 2 r, 50 

by putting the arbitrariness of the renormalization scale into the parameter 
Aqcd, known as the QCD scale parameter, In Aq C£) = In /1 2 — 16tt 2 /((]q g 2 (/i 2 )). 
Once the moments of the structure function or parton distribution are known 
at /x 2 , Eq.(|4^) can be used to give the moments at any other value of Q 2 . 

An intuitive and mathematically equivalent_fjicture for this Q 2 evolution 
is provided by the DGLAP evolution equations EJ. For the non-singlet quark 
distribution one has: 

dq NS {x,Q 2 ) a s {Q 2 ) f 1 dy 

— <lNS{y,t) Pqq(x/y) , (51) 



dt 2ir J x y 

where P qq {x/y) is the q — > q + g splitting function, which gives the probabil- 
ity of finding a quark with momentum fraction x inside a parent quark with 
momentum fraction y, after it has radiated a gluon. The splitting function is 
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related to the anomalous dimension by j N ~ J dzz N ~ 1 P(z). The generaliza- 
tion to singlet evolution is 3 E a l \ n .straightforward, but involves a set of coupled 
quark and gluon equations EfE3H. 

The physical interpretation of the DGLAP equations is that as Q 2 in- 
creases the quarks radiate more and more gluons, which subsequently split into 
quark and antiquark pairs, which themselves then radiate more gluons, and 
so on. In this manner the quark-antiquark sea can be generated from a pure 
valence component. This process modifies the population density of quarks as 
a function of x, so that the momentum carried by quarks is no longer a static 
property of the nucleon, but now depends on the resolving power of the probe, 
Q 2 . In general, the larger the Q 2 , the better the resolution, and the more 
substructure seen in the hadron. It is a remarkable success of perturbative 
QCD that it can provide a quantitative description of the scaling violations 
of structure functions 32-34 for a large range of x and over many orders of 
magnitude of Q 2 . 



3.3 Semi-Inclusive Scattering 

Inclusive electron-nucleon scattering is a well-established tool which has been 
used to study nucleon structure for many years. Somewhat less exploited, but 
potentially more powerful, is semi- inclusive scattering cB, in which a specific 
hadron, h, is observed in coincidence with the scattered electron, eiV — > e'hX. 
This process offers considerably more freedom to explore the individual quark 
content of the nucleon than is possible through inclusive scatteringl23. 

A central assumption in semi-inclusive DIS is that at high energy the 
interaction and production processes factorize. Namely, the interaction of the 
virtual photon with a parton takes place on a much shorter space-time scale 
than the fragmentation of the struck quark, and the spectator quarks, into final 
state hadrons. Furthermore, the hadronic products of the scattered quark (in 
the current fragmentation region, along the direction of the current in the 
photon-nucleon center of mass frame) should be clearly separated from the 
hadronic remnants of the target (in the target fragmentation region) . 

The cleanest way to study fragmentation is in the current fragmentation 
region, where the scattered quark fragments into hadrons by picking up qq 
pairs from the vacuum. The production of a specific final state hadron, h, is 
parameterized by a fragmentation function, D^(z), which gives the probability 
of quark q fragmenting into a hadron h with a fraction z of the quark's (or, 
at high energy, the photon's) center of mass energy. Because it requires only 
a single qq pair, the leading hadrons in this region are predominantly mesons. 
At large z, where the knocked out quark is most likely to be contained in 



19 



the produced meson, one can obtain direct information on the momentum 
distribution of the scattered quark in the target. At small z this information 
becomes diluted by additional qq pairs from the vacuum which contribute to 
secondary fragmentation. 

In the QCD-improved parton model the number of hadrons, h, produced 
at a given x, z and Q 2 can be written (in leading order) as: 

N h (x, z, Q 2 ) <*J2 e l D v( z > ■ ( 52 ) 

9 

Although factorization of the x and z dependence is generally true only at 
high energy, recent data from HERMES E3 suggests that at v ~ 10-20 GeV 
the fragmentation functionSj^re still independent of x, and agree with previous 
measurements by the EMCEj at somewhat larger energies. Where the factor- 
ization hypothesis breaks down is not known, and the proposed semi-inclusive 
program at an energy upgraded CEBAF at Jefferson Lab, with v typically ~ 
5-10 GeV, will test the limits of the parton interpretation of meson electro- 
production. 

3-4 Off-Forward Parton Distributions 

The nucleon's deep-inelastic structure functions and elastic form factors pa- 
rameterize fundamental information about its quark substructure. Both re- 
flect dynamics of the internal quark wave functions describing the same phys- 
ical ground state, albeit in different kinematic regions. An example of how 
in certain cases these are closely related is provided by the phenomenon of 
quark-hadron duality (Section 4.1). 

Recently it has been realized that form factors and structure functions 
can be simultaneously embedded within the general framework of off-forward 
(sometimes also referred to as non-forward, or skewed) parton distributionsLEEj. 
The off-forward parton distributions (OFPDs) generalize and interpolate be- 
tween the ordinary parton distributions in deep-inelastic scattering and the 
elastic form factors. 

As illustrated in Fig. 3, the OFPD is the amplitude (in the infinite momen- 
tum frame) to remove a parton with momentum k^ from a nucleon of momen- 
tum P^, and insert it back into the nucleon with momentum k'^ = fc M + P^ — P M , 
where P^ is the final state nucleon momentum. The simplest physical process 
in which the OFPD can be measured is deeply- virtual Compton scattering cB 
(DVCS). 

For the spin-averaged case, the OFPDs can be defined as matrix elements 
of bilocal operators ip(—Xn/2)jC. , y ,J 'ip(Xn/2), where A is a scalar parameter, and 
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p p' 

Figure 3: Leading twist off-forward parton distribution, as seen in deeply-virtual Compton 
scattering. The nucleon, parton and photon initial (final) momenta are labeled P (P 1 ), k 
(k r ) and q (q 1 ), respectively. Deep-inelastic scattering corresponds to q = q' , P = P' . 



n M is a light-like vector proportional to (1; 0, 0, —1). The gauge link C, which 
is along a straight line segment extending from one quark field to the other, 
makes the bilocal operator gauge invariant. In the light-like gauge A^n^ = 0, 
so that the gauge link is unity. The most general expression for the leading 
contributions at large Q 2 can then be written 

^e iXx (P'\lp{-Xn/2)Y l ^(\n/2)\P) = H(x,£,t) u(P')Y l u{P) 

+ E(x, g, t) u(P')i*^ {K 2 J u) u(P) + ■ ■ , (53) 

where t — [P 1 — P) 2 and £ = — n ■ (P' — P), with u(P) the nucleon spinor, 
and the dots (• • •) denote higher- twist contributions. A similar expression 
can be derived for the axial vector current. The structures in Eqs.(|53|) are 
identical to those in the definition of the nucleon's elastic form factors, Eq. (|l6|) . 
The chiral-even distribution, H , survives in the forward limit in which the 
nucleon helicity is conserved, while the chiral-odd distribution, E, arises from 
the nucleon helicity flip associated with a finite momentum transfer LJ. 

The off-forward parton distributions display characteristics of both the 
forward narton distributions and nucleon form factors. In the limit of P^ — > P^, 
one finds E3: 

H(x,0,0) = q(x) , (54) 



where q(x) is-the forward quark distribution, defined through similar light-cone 
correlationstil (the dependence on the scale, Q 2 , of H and q is suppressed). On 
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the other hand, the first moment of the off-forward distributions are related to 
the nucleon form factors by the following sum rules E3: 




(55) 



J-i 




(56) 



where the x-integrated distributions are in fact independent of £. These sum 
rules provide important constraints on any model calculation of the OFPDsc3. 

Higher moments of the OFPDs can also be related to matrix elements of 
the QCD energy-momentum tensor. Because the form factors of the energy- 
momentum tensor contain information about the quark and gluon contribu- 
tions to the nucleon angular momentum, the OFPDs can therefore provide 
information on the fraction of the nucleon spin carried by quarks and gluons, 
which has been a subject of intense interest now for more than a decade (see 
Section 4.4). 

Having introduced the tools necessary to study nucleon substructure in 
inclusive and exclusive reactions, in the next Section we examine more closely 
the dependence on flavor and spin of the quark momentum distributions. 

4 Flavor and Spin Content of the Nucleon 

The distribution of quarks in the nucleon is perhaps the most fundamental 
problem in hadron physics. Knowing the total structure functions and form 
factors as a function of x and Q 2 is important for determining the scaling 
properties and global characteristics of the nucleon, however, understanding 
the relative contributions from different quark flavors gives us deeper under- 
standing of the nucleon's internal structure and dynamics. In this Section we 
first explore the flavor dependence of the valence quark distributions at large 
x, then discuss the flavor structure of the quark-antiquark sea. Finally, we 
review several topical issues concerning the spin structure of the nucleon. 

4-1 Valence Quarks 

Much of the emphasis in recent years has been placed on exploring the region of 
small Bjorken-ir at high-energy colliders such as HERA. Delving into the small- 
x region is necessary in order to determine integrals of structure functions and 
minimize x — > extrapolation errors when testing various integral sum rules. 

At small x (x < 0.2) most of the strength of the structure function is due 
to the quark-antiquark sea generated through perturbative gluon radiation 
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and subsequent splitting into quark-antiquark pairs, g — > qq. Genuine non- 
perturbative effects associated with the nucleon ground state structure are 
therefore more difficult to disentangle from the perturbative background. 

Valence quark distributions, on the other hand, reflect essentially long- 
distance, or non-perturbative, aspects of nucleon structure, and can be more 
directly connected with low energy phenomenology 24-26 associated with form 
factors and the nucleon's static properties. After many years of structure 
function measurements over a range of energies and kinematical conditions, the 
valence quark structure has for some time now been thought to be understood. 
However, there is one major exception — the deep valence region, at x > 0.7. 

Knowledge of quark distributions at large x is essential for a number of 
reasons. Not least of these is the necessity of understanding backgrounds in 
collidei^xperiments, such as in searches for new physics beyond the standard 
model e3. Furthermore, the behavior of the ratio of valence d to u quark 
distributions in the limit x — ► 1 provides a critical test of the mechanism 
of spin-flavor symmetry breaking in the nucleon, and a test of the onset of 
perturbative behavior in large- a; structure functions c3. 

SU(6) Symmetry Breaking and the d/u Ratio 

The precise mechanism for the breaking of the spin-flavor SU(6) symmetry 
is a basic question for hadron structure physics. In a world of exact SU(6) 
symmetry, the wave, function of a proton, polarized say in the +z direction, 



would be given by 
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P T = A^u t (ud) s =o + t ( ud )s=i - t^u i {ud) s=1 

1 y/2 
- -d t (uu) s= i - —d i (uu) s =i , (57) 

where the subscript S denotes the total spin of the two-quark component. 
In this limit, apart from charge and isospin, the u and d quarks in the proton 
would be identical, and the nucleon and A would, for example, be degenerate in 
mass. In deep-inelastic scattering, exact SU(6) symmetry would be manifested 
in equivalent shapes for the valence quark distributions of the proton, which 
would be related simply by u va i(x) = 2d va \(x) for all x. From Eq.(|3^), for the 
neutron to proton structure function ratio this would imply: 

pn 9 

= - [SU(6) symmetry]. (58) 

In nature spin-flavor SU(6) symmetry is, of course, broken. The nucleon 
and A masses are split by some 300 MeV. Furthermore, it is known that the d 
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quark distribution in DIS is considerably softer than the u quark distribution, 
with the neutron/proton ratio deviating at large x from the SU(6) expectation. 
The correlation between the mass splitting in the 56 baryons and the large- a; 
behavior, of F^/F^ was observed some time ago 46-48 . Based on phenomeno- 
logicalhj and Reggec3 arguments, the breaking of the symmetry in Eq.(f37|) 
was argued to arise from a suppression of the 'diquark' configurations having 
5 = 1 relative to the 5 = configuration, namely: 

(qq)s=o > (qq)s=i, x -> 1 . (59) 



Such a suppression is in fact quite natural lla£3 if one observes that whatever 
mechanism leads to the observed N — A splitting (e.g. color-magnetic force, 
instanton-induced interaction, pion exchange), it necessarily acts to produce 
a mass splitting between the two possible spin states of the two quarks which 
act as spectators to the hard collision, (qq)s, with the 5 = 1 state heavier 
than the 5 = state by some 200 MeV. From Eq.(^?]), a dominant scalar 
valence diquark component of the proton suggests that in the x — ► 1 limit F^ 
is essentially given by a single quark distribution (i.e. the it), in which case: 
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> [5 = dominance] . (60) 



This expectation has, in fact, been built into most phenomenological fits to 
the parton distribution data 32-34 . 

An alternative suggestion, based on perturbative QCD, was originally for- 
mulated by Farrar and Jackson Ell. There it was argued that the exchange 
of longitudinal gluons, which are the only type permitted when the spins of 
the two quarks in (qq)s are aligned, would introduce a factor (1 — a;) 1 / 2 into 
the Compton amplitude — in comparison with the exchange of a transverse 
gluon between quarks with spins anti-aligned. In this approach the relevant 
component of the proton valence wave function at large x is that associated 
with states in which the total 'diquark' spin projection, S z , is zero: 

(qq)s z =o > (qq)s,=i, x 1 . (61) 

Consequently, scattering from a quark polarized in the opposite direction to 
the proton polarization is suppressed by a factor (1 — x) relative to the helicity- 
aligned configuration. 

This is related to the treatment based on counting rules where the large- a; 
behavior of the parton distribution for a quark polarized parallel (A5 Z — 1) or 
antiparallel (A S z = 0) to the proton helicity is given by q^ (x) = (1 — x ) 2n ~ 1+AS ' 
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where n is the minimum number of non-interacting quarks (equal to 2 for the 
valence quark distributions). In the x — > 1 limit one therefore predicts: 

F n 3 d 1 

— =■ — ► - , ► - [5 2 = dominance] . (62) 

1<2 I U 

Similar predictions can be made for the ratios of polarized quark distributions 
at large x (see Section 4.4). 

The biggest obstacle to an unambiguous determination of d/u at large x 
is the absence of free neutron targets. In practice essentially all information 
about the structure of the neutron is extracted from light nuclei, such as the 
deuteron. The deuteron cross sections must however be corrected for nuclear 
effects in the structure function, which can become quite significant E3 at large 
x. In particular, whether one correctsJjor Fermi motion only, or in addition for 
binding and nucleon off-shell effects £3, the extracted F% /F% ratio can differ 
by up to ~ 50% for beyond x ~ 0.5. 

A number of suggestions have been made how to avoid the nuclear con- 
tamination problem 53-57 . One of the more straightforward ones is to measure 
relative yields of 7r + and 7r~ mesonsjii semi-inclusive scattering from protons 
in the current fragmentation region At large z (z being the energy of the 
pion relative to the photon) the u quark fragments primarily into a 7r + , while 
a d fragments into a n~ , so that at large x and z one has a direct measure of 
d/u. 

From Eq. ( |52"| ) the number of charged pions produced from a proton target 
per interval of x and z is at leading order in QCD given byE 5 !: 

Nf (x, z, Q 2 ) ~ Au(x, Q 2 ) D{z, Q 2 ) + d(x, Q 2 ) D{z, Q 2 ) , (63) 
Nf {x, z, Q 2 - 4u(x, Q 2 ) D{z, Q 2 ) + d(x, Q 2 ) D(z, Q 2 ) , (64) 

where D = D 7 /^ = DJ is the leading fragmentation function (assuming 
isospin symmetry), and D(z) = = _DJ is the non-leading fragmenta- 

tion function. Taking the ratio of these one finds: 

V ; N* + A + d/u-D/D K ' 

In the limit z — > 1, the leading fragmentation function dominates, D D, 
and the ratio IT — > (l/4)d/u. 

In the realistic case of smaller z, the D/D term in contaminates the 
yield of fast pions originating from struck primary quarks, diluting the cross 
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section with pions produced from secondary fragmentation by picking up extra 
qq pairs from the vacuum. Nevertheless, one can estimate the yields of pions 
using the empirical fragmentation functions measured by the HERMES Col- 
laboration Ell and the EMCE3. Integrating the differential cross section over 
a range of z, as is more practical experimentally, the resulting ratios for cuts 
of z > 0.3 and z > 0.5 are shown in Fig. 4 at Q 2 ~ 5 GeV 2 for two different 
asymptotic x — > 1 behaviorscJcJ'EJ: d/u — > (dashed) and d/u — ► 1/5 (solid). 




Figure 4: Semi-inclusive pion ratio R n as a function of x for fixediz ft! 1. The dashed line 
represents the ratio constructed from the CTEQ parameterization! 32 !, while the solid includes 
the modified d distribution. 

The HERMES Collaboration has previously extracted the dhi ratio from 
the 7r + -7r - difference using both proton and deuteron targets 123 to increase 
statistics. The advantage of using both p and d is that all dependence on 
fragmentation functions cancels, removing any uncertainty that might be in- 
troduced by incomplete knowledge of the hadronization process. On the other 
hand, at large x one still must take into account the nuclear binding and Fermi 
motion effects in the deuteron, and beyond x ~ 0.7 the difference between the 
ratios with-corrected for nuclear effects and those which are not can be quite 
dramatic L3. Consequently a d/u ratio obtained from such a measurement 
without nuclear corrections could potentially give misleading results. 

On the other hand, with the high luminosity electron beam available at 
CEBAF, one will be able to compensate for the falling production rates at large 
x and z, enabling good statistics to be obtained with protons alone. Such a 
measurement will become feasible with an upgraded 12 GeV electron beam, 
which will enable greater access to the region of large Q 2 and W 2 . Since the 
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x — > 1 behavior is one of the very few predictions for the ^-dependence of 
quark distributions which can be drawn directly from QCD, the results of such 
measurements would clearly be of enormous interest. 

Quark-Hadron Duality 

Quark-hadron duality provides a beautiful illustration of the connection be- 
tween structure functions and nucleon resonance form factors 59-63 . In partic- 
ular, it allows the behavior of valence structure functions in the limi t x — ► 1 
to be determined from the Q 2 dependence of elastic form factors 

The subject of quark-hadron duality, and the relation between exclusive 
and inclusive processes, is actually as old as the first deep-inelastic scattering 
experiments themselves. In the early 1970s the inclusive-exclusive connection 
was studied in the context of deep-inelastic scattering in the resonance region 
and the—onset of scaling behavior. In their pioneering analysis, Bloom and 
GilmanEj observed that the inclusive F2 structure function at low W gener- 
ally follows a global scaling curve which describes high W data, to which the 
resonance structure function averages. Furthermore, the equivalence of the 
averaged resonance and scaling structure functions appears to hold for each 
resonance, over restricted regions in W, so that the resonance — scaling duality 
also exists locallyEJ. 

Following-JJloom and Gilman's empirical observations, de Rujula, Georgi 
and PolitzerEil pointed out that global duality can be understood from an op- 
erator product expansion of QCD moments of structure functions. Expanding 
the F2 moments in a power series in 1/Q 2 (c.f. Eq.(|39j)), 

f Q d£ i N - 2 F 2 (g, Q 2 ) = g ( {N Q2 2) ^ ) " Af (a s (Q 2 )) , (66) 

where A is some mass scale, and the Nachtmann scaling variable £ = 2a:/ (1 + 
yl + x 1 fr) takes into account target mass corrections, one can attribute the 
existence of global duality to the relative size of higher twists in deep-inelastic 
scattering. The Q 2 dependence of the coefficients Ajy arises only through 
a s (Q 2 ) corrections, and the higher twist matrix elements A^ >0 ' are expected 
to be of the same order of magnitude as the leading twist term A$ . The weak 
Q 2 dependence of the low F2 moments can then be interpreted as indicating 
that higher twist (l/Q 2k suppressed) contributions are either small or cancel. 

Although global Bloom-Gilman duality of low structure function moments 
can be analyzed systematically within a perturbative operator product expan- 
sion, an elementary understanding of local duality's origins is more elusive. 
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This problem is closely related to the question of how to build up a_scaling 
(« Q 2 independent) structure function from resonance contributions E3, each 
of which is described by a form factor G R (Q 2 ) that falls off as some power of 



To illustrate the interplay between resonances and scaling functions, one 
can observe EjI£3 that (in the narrow resonance approximation) if the contri- 
bution of a resonance of mass Mr to the F2 structure function at large Q 2 is 



given by (c.f. Eqs.@-@) F 2 (ii) = 2Mv (G R (Q 2 )) 2 S(W 2 - M|), then a 



form factor behavior G R (Q 2 ) ~ (l/Q 2 )™ translates into a structure function 



F^ H) - (1 - x R ) 2n ~\ where x R = Q 2 /(M 2 - M 2 + Q 2 ). On purely kine- 



matical grounds, therefore, the resonance peak at x R does not disappear with 
increasing Q 2 , but rather moves towards x = 1. 

For elastic scattering, the connection between the 1/Q 2 power of the elas- 
tic form factors at large Q 2 and the. a; — > 1 behavior of structure functions was 
first established by Drell and YanEil and WestEHl. Although derived before the 
advent of QCD, the Drell- Yan — West form factor-structure function relation 
can be expressed in perturbative QCD language in terms of hard gluon ex- 
change. The pertinent observation is that deep-inelastic scattering at x ~ 1 
probes a highly asymmetric configuration in the nucleon in which one of the 
quarks goes far off-shell after exchange of at least two hard gluons in the initial 
state; elastic scattering, on the other hand, requires at least two gluonsJxi the 
final state to redistribute the large Q 2 absorbed by the recoiling quark t3. 

If the inclusive-exclusive connection via local duality is taken seriously, 
one can use measured structure functipns in the resonance region at large £ to 
directly extract elastic form factors E3. Conversely, empirical electromagnetic 
form factors at large Q 2 canJpe used to predict the x — > 1 behavior of deep- 



Integrating the elastic contributions to the structure functions in Eqs.(|19|)- 
( p2| ) over the Nachtmann variable £, where £ = 2x/(l + y/l + x 2 /t), between 
the pion threshold £th and £ = 1, one finds 'localized' moments of the structure 
functions: 



I/O 2 . 





(67) 



(68) 




N-2 



4 - 2£ 1 + t 



^ Gm{Q 2 )(G e {Q 2 )+tG m {Q 2 )) 



, (69) 
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/■ « - -4- ^wi(y-c.OT} , (70) 

4 h 4 - 2£ 1 + r 

where r = Q 2 /4M 2 and Co = 2/(1 + yj\ + 1/r) is the value of f at x = 1. 
Differentiating Eqs. j67|)— (|70|) with respect to Q 2 for N = 2 allows the inclusive 
structure functions near^c = 1 to be extracted from the elastic form factors 
and their Q 2 -derivativesL3: 

« « ^ . Cn> 

Gm{Gm — Ge) , 1 ((1{GeGm) . dG M \ 

gl0C 4jif»(i + r)» + TT7( v ^g^ + T dg^J ' (73) 

Gm(Gm — Ge) , t (<1(GeGm) dG 2 j\ 
92(X AM 2 {l + r) 2 + —{ dQ 2 • (?4) 

Note that as r — > oo each of the structure functions Pi , P2 and gi is determined 
by the slope of the square of the magnetic form factor, while 52 (which in 
deep-inelastic scattering is associated with higher twists) is determined by a 
combination of Ge and Gm- 

Equations (|7l"|)-(|7"4|) allow the x ~ 1 behavior of structure functions to 
be predicted from empirical electromagnetic form factors. The ratios of the 
neutron to proton P\, P 2 and g\ structure functions are shown in Fig. 5 as 
a function of Q 2 , using typical parameterizations EJ of the global form fac- 
tor data. While the P2 ratio varies quite rapidly at low Q 2 , beyond Q 2 ~ 
3 GeV 2 it remains almost Q 2 independent, approaching the asymptotic value 
(dG\ 2 1 'dQ 2 ) I \dG M I 'dQ 2 ) . This is consistent with the operator product expan- 
sion interpretation of de Riijula et al. in which duality should be a better 
approximation with increasing Q 2 . Because the F™/Ff ratio depends only on 
Gm, it remains flat over nearly the entire range of Q 2 . At asymptotic Q 2 the 
model predictions for F\{x — > 1) coincide with those for F2; at finite Q 2 the 
difference between Pi and P2 can be used to predict the x — ► 1 behavior of the 
longitudinal structure function, or the R = <Jl/&t ratio. 

The pattern of SU(6) breaking for the spin-dependent structure function 
ratio g\j g\ essentially follows that for F^/Ft?, namely l/4in_trie d quark 
suppression and 3/7 in the helicity flip suppression scenarios EjH However, 
the 51 structure function ratio approaches the asymptotic limit somewhat more 
slowly than Pi or P2, which may indicate a more important role played by 
higher twists in spin-dependent structure functions than in spin-averaged. 
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It appears to be an interesting coincidence that the helicity retention model 
prediction of 3/7 is very close to the empirical ratio of the squares of the neutron 
and proton magnetic form factors, ~ 4/9. Indeed, if one approximates 

the Q 2 dependence of the proton and neutron form factors by dipoles, and 
takes G E ~ 0, then the structure function ratios are all given by simple analytic 
expressions, F^/F^ « F^/Ff w ffi/ffi ~ > Mn/^p as Q 2 — * °°- On the other 
hand, for the 52 structure function, which depends on both Ge and Gm at large 
Q 2 , one has a different asymptotic behavior, g% j g\ — > /i^/(/i p (l + /i p )) w 0.345. 
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Figure 5: Neutron to proton ratio for F\ (dashed), F2 (solid) and gi (dot-dashed) structure 
functions in the limit x — > 1. 

If the resonance structure functions at large £ are known, one can con- 
versely extract the nucleon electromagnetic form factors from Eqs.(|6^)-([70|). 
The Gm form factor of the nucleon can be extracted directly from the measured 
-Fi(£iQ 2 ) structure function in Eq.(|67|). Unfortunately, only the i*2(£, Q 2 ) 
structure function of the proton has so far been measured in the resonance 
region. Nevertheless, to a good approximation one can assume that the ra- 
tio of-filectric to magnetic form factors is reasonably well known (see, however, 
Ref.Ea), and extract Gm from the F2 structure function in the resonance region 
via Eq.(ff). 

Using the parameterization of the recent ^(CQ 2 ) data from Jefferson 
LabE3, in Fig. 6 we show the extracted G M compared with a compilation of 
elastic data. The agreement with data is quite remarkable over the entire range 
of Q 2 between and 3 GeV 2 . 

The reliability of the duality predictions is of course only as good as the 
quality of the empirical data on the electromagnetic form factors and reso- 
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Figure 6: Proton magnetic form factor extracted from the inclusive structure function via 
Eq.@. 

nance structure functions. While the duality relations Efl are expected to be 
progressively more accurate with increasing Q 2 , the difficulty in measuring 
form factors at large Q 2 also increases. Experimentally, the proton magnetic 
form factor G V M is relatively well constrained to Q 2 ~ 30 GeV 2 , and the pro- 
ton electric G P E to Q 2 ~ 10 GeV 2 . The neutron magnetic form factor 
has been measured to Q 2 ~ 5 GeV 2 , although the neutron G% is not very 
well determined at large Q 2 (fortunately, however, this plays only a minor role 
in the duality relations, with the exception of the neutron to proton gi ratio, 
Eq.@). 

Obviously more data at larger Q 2 would allow more accurate predictions 
for the x — > 1 structure functions, and new experiments at Jefferson Lab H 
will provide valuable constraints. Once data on the longitudinal and spin- 
dependent structure functions at large x become available, a more complete 
test of local duality between elastic form factors and x ~ 1 structure functions 
can be made. 

4-2 Light Quark Sea 

Over the past decade a number of high-energy experiments and refined data 
analyses have forced a re-evaluation of our view of the nucleon in terms of 
three valence quarks immersed in a sea of perturbatively generated qq pairs 
and gluonsEa. A classic example of this is the asymmetry of the light quark sea 
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of the proton, dramatically confirmed in._Eecent deep-inelastic and Drell-Yan 
experiments at CERNEJO and FermilabO. 

Difference between quark or antiquark distributions in the proton sea al- 
most universally signal the presence of phenomena which require understand- 
ing of strongly coupled QCD. Their existence testifies to the relevance of long- 
distance dynamics (which are responsible for confinement) even at large energy 
and momentum transfers. 

Because gluons in QCD are flavor-blind, the perturbative process g — > qq 
gives rise to a sea component of the nucleon which is symmetric in the quark 
flavors. Although differences can arise due to different quark masses, because 
isospin symmetry is such a good symmetry in nature, one would expect that 
the sea of light quarks generated perturbatively would be almost identical, 
u{x) = d{x). 

It was therefore a surprise to_many when measurements by the New Muon 
Collaboration (NMC) at CERNLJ of the proton and deuteron structure func- 
tions suggested a significant excess of J over u in the proton. Indeed, it heralded 
a renewed interest in the application of ideas from non-perturbative QCD to 
deep-inelastic scattering analyses. While the NMC experiment measured the 
integral of the antiquark difference, more recently the E866 experiment at Fer- 
milab has for the first time mapped out the shape of the d/u ratio over a large 
range of x, 0.02 < x < 0.345. 

Specifically, the E866/NuSea Collaboration measured Drell-Yan pairs 

produced in pp and pd collisions. In the parton model the Drell-Yan cross sec- 
tion is proportional to: 

o» h oc £ <^) + ?(xi) q h (x 2 )) , (75) 

where h = p or D, and x\ and xi are the light-cone momentum fractions 
carried by partons in the projectile and target hadron, respectively. Using 
isospin symmetry to relate quark distributions in the neutron to those in the 
proton, in the limit x\ 3> x 2 ( m which q{x\) <C q{x\f) the ratio of the deuteron 
to proton cross sections can be written: 

a* _ 1 ( 1 d(x 2 ) \ A + d( Xl )/u( Xl ) 

2<jpp 2\ u{x 2 )) 4 + d(a?i)/u(x 1 ) • d(x 2 )/u(x 2 ) ' 1 ' 

Corrections for nuclear shadowing in the deuteron El-El, which are important 
at x <C 0.1, are small in the region covered by this experiment. 

The relatively large asymmetry found in these experiments, shown in Fig. 7 
implies the presence of non-trivial dynamics in the proton sea which does not 
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have a perturbative QCD origin. The simplest and most obvious source of a 
non-perturbative asymmetry in the light quark sea is the chiral structure of 
QCD. From numerous studies in low energy physics, including chiral pertur- 
bation theory E£l, pions are known to play a crucial role in the structure and 
dynamics of the nucleon (see Section 2). However, there is no reason why the 
long-range tail of the nucleon should not also play a role at higher energies. 




Figure 7: Flavor asymmetry of the light antiquark sea, including pion cloud (dashed) and 
Pauli blocking effects (dotted), and the total (solid). 

As pointed out by Thomas 0, if the proton's wave function contains an 
explicit 7r + n Fock state component, a deep-inelastic probe scattering from the 
virtual 7r + , which contains a valence d quark, will automatically lead to a d 
excess in the proton. To be specific, consider a model in which the nucleon 
core consists of valence quarks, interacting via gluon exchange for example, 
with sea quark effects introduced through the coupling of the core to qq states 
with pseudoscalar meson quantum numbers (many variants of such a model 
exist — see for example RefsL3'LZl) . The physical nucleon state (momentum P) 
can then be expanded (in the one-meson approximation) as a series involving 
bare nucleon and two-particle meson-baryon states: 

\N(P)) phys = Vz {|A(P)) barc 

+ I dy d 2 k ± g MN B 0BAi{y,k ± ) \B(y,k ± );M(l-y,-k x )) } , (77) 

B,M 

where M — n, K, ■ ■ ■ and B = N, A, A, • • •. The function 4>BAi(y, kj_) is the 
probability amplitude for the physical nucleon N to be in a state consisting of 
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a baryon B and meson M, having transverse momenta kj_ and — k^, and car- 
rying longitudinal momentum fractions y = k + /P + and 1— y = (P+ — k + )/P+, 
respectively. The bare nucleon probability is denoted by Z, and qmnb is the 
MNB coupling constant. The one- meson approximation in Eq.(|77j) will be 
valid as long as the meson cloud is relatively soft (Z w 1). It will progressively 
break down for harder MNB vertex functions, at which point one will need 
to include two- meson and higher order Fock state components in Eq. (|77j) . 

Because of their small masses, the irN component of the Fock state expan- 
sion in Eq.(|77j) will be the most important (contributions from heavier mesons 
and baryons will be progressively suppressed with increasing mass). In the 
impulse approximation, deep-inelastic scattering from the ttN component of 
the proton can then be understood in the IMF as the probability for a pion to 
be emitted Jay the proton, folded with the probability of finding the_a parton 
in the pionlla. For the antiquark asymmetry, this can be written asE3: 

d( x ) - u(x) = \( ^-fMv) tfix/y) > (78) 

where q* is the (valence) quark distribution in the pion (e.g. u in 7r + , normal- 
ized to unity), and the distribution of pions with a recoiling nucleon (N — > Tin 
splitting function) is given by 75 > 80 ~ 82 ; 

f UA 3, 9 ^ v f d 2 k T J* N (s« N ) f k 2 +y 2 M 2 \ 

fMv) = J {1 _ y)y (M2 _ s „ n) 2 { — y ) ■ (79) 

For the irNN vertex a pseudoscalar 175 interaction has been used, although 
in the IMF the same results are obtained with a pseudovector coupling. The 
invariant mass squared of the ttN system is given by: 

_ k\ + ml k 2 T +M 2 

SttJV — 1 : , (»Uj 

y 1 -y 

and for the functional form of the ttNN vertex function J-"7rJv(s7rAr) we can take 
a simple dipole parameterization, J~tz n (^7r n*) — ((A 2 + M 2 )/(A 2 + s^jy)) , 
normalized so that the coupling constant g^NN has its standard value (— 
13.07) at the pole (T(M 2 ) — 1). Note that the contribution from the pion 



cloud in Eq.(78) is a leading twist effect, which scales with Q 2 (at leading order 
the Q 2 dependence of d — u in Eq. ((78]) enters through the leading-twist quark 
distribution in the pion, q* . 

Ope can easily generalize the above to include higher Fock state compo- 
nents Ea, most notably the A. Because the dominant process there is p — > 
A ++ 7T~, the A will actually cancel some of the d excess generated through the 
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irN component, although this will be somewhat smaller due to the larger mass 
of the A. 

The relative contributions are partly determined by the ttNN and nNA 
vertex form factor. The form factor cut-offs A can be determined phenomenn- 
logically by comparing against various inclusive and semi-inclusive data E3, 
although the most direct way to fix these parameters is through a comparison 
of the axial form factors for the nucleon and for the N-A transition. Within 
the framework of PCAC these form factors are directly related to the corre- 
sponding form factors for pion emission or absorption. The data on the axial 
form factor are best fit, in a dipole parameterization__by a 1.3 (1-02) GeV 
dipole for the axial N (N-A transition) form factor E3, which gives a pion 
probability in the proton of w 13%(10%). 

With these parameters Fig .6 shows the d/u ratio in the proton due tp-jrA^ 
and it A components of the nucleon wave function (dashed line)E3. DataE3 on 
the sum of the u and d (which is dominated by perturbative contributions) has 
been used to convert the calculated d—u difference to the d/u ratio. The results 
suggest that with pions alone one can account for about half of the observed 
asymmetry, leaving room for possible contributions from other mechanisms. 

Another mechanism which could also contribute to the d — u asymmetry 
is associated- with the effects of antisymmetrization of qq pairs created inside 
the core E£fE3. As pointed out originally by Field and Feynman e3, because 
the valence quark flavors are unequally represented in the proton, the Pauli 
exclusion principle will affect the likelihood with which qq pairs can be created 
in different flavor channels. Since the proton contains 2 valence u quarks 
compared with only one valence d quark, uu pair creation will be suppressed 
relative to dd creation. In the ground state of the proton the suppression will 
be in the ratio d : u — 5 : 4. 

Phcnomcnological analyses in terms of low energy models (specifically, the 
MIT bag model E3) suggest that the contribution from Pauli blocking can be 
parameterized as (d — u) Pauh = r Pauh (a + 1)(1 — x) a , where a is some 
large power, with normalization, r p , less than ss 25%. Phenomenologically, 
one finds a good fit with a « 14 and a normalization r Pauh « 7%, which 
is at the lower end of the expected scale but consistent with the bag model 
predictions c3. Together with the integrated asymmetry from pions, r 71 ' ~ 0.05, 
the combined value r = T 71 " + r Pauh w 0.12 is in quite reasonable agreement 
with the experimental result, 0.100 ± 0.018 from E866. 

Although the combined pion cloud and Pauli blocking mechanisms are able 
to fit the E866 data reasonable well at small and intermediate x (x < 0.2), it 
is difficult to reproduce the apparent trend in the data at large x towards zero 
asymmetry, and possibly even an excess of u for x > 0.3. Unfortunately, the 



35 



error bars are quite large beyond x ~ 0.25, and it is not clear whether any new 
Drell-Yan data will be forthcoming in the near future to clarify this. 

A solution might be available, however, through semi-inclusive scattering, 
tagging charged pions produced off protons and neutrons. Taking the^atio of 
the isovector combination of cross sections for ir + and tt~ production^ 

iV; ++ "~ - Nf+*~ _ 3 f u-d-d + u \ ( D + D 
N% + "*~ -i\r£ +_7r ~ ~ 5 \u- d + d- u) \D-D 

the difference d — u can be directly measured provided the u and d quark 
distributions and fragmentation functions are known. The HERMES Collab- 
oration has in fact recently measured this ratioEa, although there the rapidly 
falling cross sections at large x make measurements beyond x ~ 0.3 challeng- 
ing. On the other hand, a high luminosity electron beam such as that available 
at Jefferson Lab, could, with higher energy, allow the asymmetry to be mea- 
sured well beyond x ~ 0.3 with relatively small errors. This would parallel 
the semi- inclusive measurement of the d/u ratio through Eq.(|35l) at somewhat 
larger x. 
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Figure 8: Neutron electric form factor in the pion cloud model|£3- The direct tt coupling 
contribution is labeled "tt~ p", and the recoil proton u pn~ " . 



If a pseudoscalar cloud of qq states plays an important role in the d/u 
asymmetry, its effects should also be visibleJri other flavor-sensitive observ- 
ables, such as electromagnetic form-fjactors E3. An excellent example is the 
electric form factor of the neutron Ell, a non-zero value for which can arise 
from a pion cloud, n — * pn~ . Although in practice other effects 92-94 such 
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as spin-dependent interactions due to one gluon exchange between quarks in 
the core will certainly contribute at some level, it is important nevertheless to 
test the consistency ol the above model by evaluating its consequences for all 
observables that may carry its signature. 

To illustrate the sole effect of the pion cloud, all residual interactions be- 
tween quarks in the core can be switched off, so that the form factors have 
only two contributions: one in which the photon couples to the virtual it and 
one where the photon couples to the recoil nucleon: 

FiAQ 2 ) = £dy (f$(y,Q 2 ) + f[ N 2 \y,Q 2 )) ■ (82) 
The recoil nucleon contribution is described by the functions: 

%91nN f d2k -L F{sN-n,i) F{sNirj) 



tt">(y,Q z ) 



16^ 3 J y 2 (l-y)(s N ^i-M 2 )(s N ^f-M 2 ) 

-,2 



x(ki + M'(i-vy-(i- y y2±) , (83) 



c( N )/„, rfi\ - S 9kNN f d k -L F{sN*,i) FjsNirj) 

y 2 (l - y) (s Nv>i - M 2 )(s NnJ - A/ 2 ) 
x(-2M 2 )(l-y) 2 , (84) 



where the squared center of mass energies are: 



s Niv,i(f) = SnN + — ■ [{I - y)— ± k ± J , (85) 



with Sttat defined in Eq. (pCf) . 

The contribution from coupling directly to the pion is: 



ti\y,Q 2 ) 



16tt 3 J y 2 (l-y)(s vN , l -M 2 )(s 7rN j-M 2 ) 
x(fci + M 2 (l- y ) 2 -y 2 4) > (86) 



16tt3 J y 2 {l-y){ S7! N,-M 2 ){ s „ NJ -M 2 ) 
x (2M 2 y(l-y)) , (87) 

where the irN squared center of mass energies are: 

<u ( qx , , 
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The N — > ttN splitting functions are related to the distribution functions in 
Eqs.@ and @ by: 



fr(y,Q 2 = Q) = U N (y) , 

f[ N \y,Q 2 = 0) = f Nv {y) = U N (l-y). 



(89) 
(90) 



Using the same pion cloud parameters as in the calculation of the d/u asym- 
metry in Fig. 7, the relative contributions to from the 7r _ and recoil proton 
are shown in Fig. 8. Both are large in magnitude but opposite in sign, so that 
the combined effects cancel to give a small positive G%, consistent with the 
data. Note, however, that the Pauli blocking effect plays no role in form fac- 
tors, since any suppression of u relative to d here would be accompanied by an 
equal and opposite suppression of it sea relative to d sca , and form factors always 
contain charge conjugation odd (valence) combinations of flavors. 

The fact that the model prediction underestimates the strength of the 
observed G% suggests that other mechanisms, such as the color hyperfine .in,- 
teraction generated by one-gluon exchange between quarks in the core£3£3, 
are likely to be responsible for some of the diffcrcncc-JThc lowest order Hamil- 
tonian for the color-magnetic hyperfine interaction E3 between two quarks is 
proportional to (a s /minij)Si ■ Sj. Because this force is repulsive if the spins 
of the quarks are parallel and attractive if they're antiparallel, from the SU(6) 
wave function in Eq.(j5?]) it naturally leads to an increase in the mass of the 
A and a lowering of the mass of the nucleon. The same force also leads to the 
softeningEJ of the d quark distribution relative to the u (see Eq.([59|)). Further- 
more, it leads to a distortion of the spatial (and hence charge) distributions of 
quarks in the neutron, pushing the two (negatively charged) d to the periphery 
of the neutron, while forcing the. (positively charged) u in the center, giving 
rise to a negative charge radius Ej, e » r f )„■ 

In the harmonic oscillator model at leading orcka: in a s , the hyperfine 



where a can be related to decay amplitudes and charge radiiE3. Taking the 
value a = 0.243 from the ratio of neutron to proton charge radii, the resulting 
form factor in Fig. 9 agress quite well with the available G% data. More accu- 
rate data, which will soon be available from Jefferson Lab and elsewhere over a 
range of Q 2 will allow more systematic comparison of the various mechanisms 
which contribute to SU(6) symmetry breaking. 




(91) 
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Figure 9: Electric neutron form factor in the hyperfine perturbed quark model with a har- 
monic oscillator potential. 

4-3 Strange Quarks in the Nucleon 

A complication in studying the light quark sea is the fact that non-perturbative 
features associated with u and d quarks are intrinsically correlated with the 
valence core of the proton, so that effects of qq pairs can be difficult to dis- 
tinguish from those of antisymmetrization or residual interactions of quarks in 
the core. The strange sector, on the other hand, where antisymmetrization be- 
tween sea and valence quarks plays no role, is therefore more likely to pcavide 
direct information about the non-perturbative origin of the nucleon seacS 

Evidence for non-perturbative strangeness is currently being sought in 
a number of processes, ranging from semi-inclusive neutrino induced deep- 
inelastic scattering to parity violating electron-proton scattering. As for the 
d — u asymmetry, perturbative QCD alone generates identical s and s distri- 
butions, so that any asymmetry would have to be non-perturbative in origin. 

In deep-inelastic scattering, the CCFR collaboration t3 analyzed charm 
production cross sections in v and v reactions, which probe the s and s distri- 
butions in the nucleon, respectively. The resulting difference s — s, indicated in 
Fig. 10 by the shaded area, has been extracted from the s/s ratio and absolute 
values of s + s from global data parameterizations. 

The curve in Fig. 10 corresponds to the chiral cloud model prediction for 
the asymmetry (in analogy with the pion cloud in Section 4.2), in which the 
strangeness in the nucleon is carried by kapns and hyperons, so that the s 
and s quarks have quite different origins E3E3. Taking the A hyperon as an 
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Figure 10: Strange quark asymmetry in the proton arising from a kaon cloud of the nuclieem. 
The shaded region indicates current experimental limits from the CCFR Col 1 aboratiorj P 6 I. 



illustration (the results generalize straightforwardly to other hyperons such as 
the E), the difference between the s and s can be writtenEj: 



s(x) - s(x) = I $L (f AK (y) s A (x/y) - f KA (y) s K (x/y)) , 

j x y 



(92) 



where the K distribution function jxh is the analog of the ttN splitting func- 
tion in Eq. (JZOj) , and the corresponding A distribution /ak(?/) = Ika(^ ~~ v)- 

In the IMF parameterization of the KNY (Y — A, E) vertex function, 
because the s distribution in a kaon is much harder than the s distribution 
in a hyperon, the resulting s — s difference is negative at large x, despite 
the kaon distribution in the nucleon being slightly softer than the hyperon 
distributionEa. With a dipole cut-off mass of A ~ 1 GeV, the kaon probability 
in the nucleon is « 3%. On the other hand, the exact shape and even sign 
of the s — s difference as a function of x is quite sensitive to the shape of the 
KNY vertexS. 

Overall, while the current experimental s — s difference is consistent with 
zero, it does also consistent with a small amount of non-perturbative strangeness, 
which would be generated from a kaon cloud around the nucleon E3. Of course 
other, heavier strange mesons and hyperorrs can be added to the analysis, al- 
though in the context of chiral symmetry Lj the justification for inclusion of 
heavier pseudoscalar as well as vector mesons is less clear. The addition of the 
towers of heavier mesons and baryons has also been shown in a quark modelEU 
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to lead to significant cancellations, leaving the net strangeness in the nucleon 
quite small. 




Figure 11: Strange electromagnetic form factors of the proton compared with a kaon cloud 
prediction, witlt-Hiie magnetic (M), electric (E) and total contributions indicated. For the 
HAPPEX datafefl r ss 0.4. 

Within the same formalism as used to discuss strange quark distributions 
one can also calculate the strangeness form factors of the nucleon, which are 
being .measured in parity-violating electron scattering experiments at MIT- 



Bates \) 0( \ and Jefferson Lab El The HAPPEX Collaboration at Jefferson 



LabEil has recently measured the left-right asymmetry Alr in ep — > ep elastic 



scattering, which measures the j*Z interference term: 
^ _ vr- &l ( —G F \ I 

x 



°R + <?L \ira em V2j e Gf + r G p * 



e G P E G p ^ z) + r G P M Glf ] - \{1 - 4 sin 2 6 W ) e' G P M G p } z) j , (93) 

with e = (l + 2(l + r )tan 2 ((9/2))~ 1 , and s' = ^Jt{1 + t)(\ - e 2 ) (the Q 2 
dependence in all form factors is implicit). Using isospin symmetry, one can 
relate the electric and magnetic form factors for photon and Z-boson exchange 
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where G E ^ 1S the isovector form factor. At forward angles, the asymmetry is 
sensitive to a combination of strange electric and magnetic form factors shown 
in Fig. 11 for Q 2 = 0.48 GeV 2 . With a soft KNY forpxfactor the contributions 
to both G S E and G S M are small and slightly positive in agreement with the 
trend of the data. This result is consistent with the earlier experiment by the 



SAMPLE Collaboration at MIT-Bates 1123, G S M = +0.23 ± 0.44 at Q 2 = 0.1 



GeV 2 in a similar experiment but at backward angles. 

Although the experimental results on non-perturbative strangeness in both 
structure functions and form factors are still consistent with zero, they are 
nevertheless compatible with a soft kaon cloud around the nucleon. Future 
data on G S E , M from the HAPPEX-II and GO experiments at Jefferson Lab with 

smaller error bars and over a large range of Q 2 , as well as the remaining data 
on the proton and deuteron from SAMPLE, will hopefully provide conclusive 
evidence for the presence or otherwise of a tangible non-perturbative strange 
component in the nucleon. 



4-4 Polarized Quarks 

Most of the discussion thus far has dealt with the flavor dependence of quarks 
in the nucleon. On the other hand, there has been considerable interest over 
the past decade in how the spin of the nucleon is distributed amongst its con- 
stituents 102-104 . Spin degrees of freedom allow access to information about the 
structure and interactions of hadrons which would otherwise not be available 
through unpolarizcd processes. Indeed, experiments involving spin-polarized 
beams and targets have often yielded surprising results and presented severe 
challenges to existing theories. 

A fundamental sum rule for the spin in the nucleon states that: 

~ = Jg + Jg > (95) 

where J q and J g are the total quark and gluon angular momenta, which can 
be decomposed into their helicity and orbital contributions: 

J q = ^AE + L q , (96) 

Jg = AG + Lg . (97) 

In particular, AS, which is defined as the forward matrix element of the ax- 
ial current, AE = (N\'tpj^5^>\N), measures the total helicity of the nucleon 
carried by quarks, which for three flavors is: 

AT, = Au + Ad + As . (98) 
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In non-relativistic quark models the spin of the nucleon is carried entirely by 
valence quarks, so that ASnrqm = 1< 

The gluon helicity, AG, can be measured in high-energy polarized proton- 
proton collisions, via charm production through 



ark-gluon fusion, or in pro- 
duction of jets with high transverse momentum tH3. The orbital contributions, 
L q and L g , can in principle be extracted from measurements of off- forward pax- 
ton distributions in DVCS, or deeply-virtual meson production experimentsLJ. 
Currently there is little empirical information on the gluon polarization, and 
the angular momentum distributions are totally unknown. Note that each 
term in Eqs.(^6j) and (|97]) is renormalization scheme and scale dependent, and 
only the quark helicity can be defined in a gauge- invariant manner tj. 

Experimentally, AE (which is also referred to as the singlet axial charge) 
can be determined from a combination of triplet and octet axial charges, 
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which are determined from /3-decays of the nucleon and hyperons, and the 
spin-dependent g\ structure function of the nucleon. In the MS scheme, the 
lowest moment of-the-jC^— ptructure function (at lowest order in perturbative 
QCD) is given byt 



104, LOf . 



dxg{ (n) (at, 



1 



(t) 



4* + S» + 5 AE ) • (101) 
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where the ± refers to the proton or neutron. 

The first spin structure function experiments at CERN 
rather small value for AS, in fact consistent with zero, which prompted the 
so-called 'proton spin-crisis'. A decade of subsequent measurements of inclu 
sive spin structure functions using proton, deuteron and 3 He targets have 
determined AT, much more accurately, with the current world average value 
(at a scale of Q 2 = 10 GeV 2 in the MS scheme) being AE w 0.3. 

While the spin fractions carried by quarks in the nucleon require only the 
first moment of the inclusive spin-dependent structure functions, to determine 
the x dependence of the polarized distributions requires independent linear 
combinations of Aq which at present can only be obtained from semi-inclusive 
scattering. Generalizing Eq.(^2|) to the production of hadrons with a polarized 
beam and target, eN — > e'hX, the difference between the number of hadrons 
produced at a given x, z and Q 2 with electron and nucleon spins parallel and 
antiparallel is (at leading order) given by: 



AN h = N^_ lf (x,z,i 



2 ) 



Aq(x,Q 2 ) AD h Jz,Q 2 



(102) 
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where the polarized fragmentation function AD^ gives the probability for a 
polarized quark q to hadronize into a hadron h. 

In analogy with the large- a; behavior of the unpolarized u and d distribu- 
tions in Section 4.1, the x — > 1 limit of polarized quark distributions provides 
a sensitive test of various mechanisms of spin-flavor symmetry breaking. For 
SU(6) symmetry, the ratio of the polarized to unpolarized quark distributions 
is: 

^ = I , ^ = ~\ [SU(6) symmetry] . (103) 

If the symmetry is broken through the suppression of the S = 1 diquark 
contributions in the nucleon, then in the limit x — > 1: 

— — > 1 , — ^ — ► -- [5 = dominance] . (104) 
u a 3 

The perturbative QCD prediction (where the dominant configurations of the 
proton wave function are those in which the spins of the interacting quark and 
proton are aligned) on the other hand, is: 

— — ► 1 , ^ — > 1 [S z = dominance] . (105) 
u a 

Note that the predictions for the Ad quark in particular are quite different in 
the perturbative and non-perturbative models, even differing by sign. 

The spin-flavor distributions can be directly measured via polarization 
asymmetries fpt-the. difference between tt + and 7r~ production cross sections 
on the protonL 1 



lOS.LlO 



... AN; - AN; 4A Mval - Ad val 

AN£ + +AN£~ 4u va i-d va i ' ' 



where the dependence on fragmentation functions and sea quarks cancels. A 
combination of inclusive and semi-inclusive asymmetries using protons and 
deuterons at an energy upgraded Jefferson Lab will allow the spin-dependent 
Au and Ad distributions to be determined up to x ~ 0.8 with good statis- 
tics liiil, which sho ul d be able to discriminate between the various model sce- 
narios in Eqs.(fa|-(|l0f). 



At smaller x, similar combinations of asymmetries could also be able to 
measure the polarized antiquark distributions, Ad and Ail. These are par- 
ticularly interesting in view of the qualitatively different predictions in non- 
perturbative models. While chiral (pion) cloud models do not allow any po- 
larization in the antiquark sea, the Pauli exclusion principle on the other hand 
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predicts quite a large asymmetry, (Au — Ad) / (d— u) — 5/3, even bigger than in 
the unpolarized seacJiH3. Measurement of the polarized asymmetry in semi- 
inclusive scattering would then enable the relative sizes of the pion and Pauli 
blocking contributions to d — u to be disentangled. 

While data on the polarized quark distributions has slowly been accumu- 
lating from various experiments, and plans are under way to systematically 
measure the polarization of the gluons, until recently there has been very little 
discussion about the fraction of the nucleon spin residing in_angular momen- 
tumEHl. This changed somewhat when it was demonstrated!^ that the orbital 



angular momentum contributions could be determined from off-forward parton 
distributions measured in deeply-.sd.rtual Compton scattering (see Section 3.4). 

In particular, it was shown 123 that a sum rule can be derived relating 
moments of the OFPDs to the total angular momentum carried by quarks and 
gluons: 



dx x(H(x, f , t) + E(x, £, t)) = A(t) + B(t) , (107) 

where A and B are form factors of the energy-momentum tensor in QCD0: 

T *" = ~ (j J ^ ir ^^ + ^l^iD^^) , (108) 

(109) 



9 ~ 4 r a ■ 

for quarks and gluons, respectively, where the braces {• • •} represent_sym- 
metrization of indices. The matrix elements of can be expanded asEj: 



(P'\T" V \P) = u(P') J-A(t) 7 { ^(P + P'P 

+ lM B(t) (P + P ') { ^ }a ( P '- p )c 



u(P) . (110) 



One can then show that the total angular momentum carried by quarks is 
given by: 



J q = l(A(Q)+B(0) 



(111) 



Combining the extracted J q with AS measured in inclusive DIS, one can then 
determine the orbital angular momentum of the quarks in the nucleon. An 
analogous sum rule can also be written for the total gluon angular momentum, 
J g , which can be obtained from OFPDs measured in deeply- virtual meson 
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production. From this the gluon orbital angular momentum can be extracted 
once the gluon helicity AG is known. 

The program to measure the off-forward parton distributions H,E, - ■ ■ in 
deeply-virtual Compton scattering and meson production experiments is dif- 
ficult, requiring a large coveragej-of. kinematics and knowledge of background 



such as the Bethe-Heitler processEj for DVCS. The first steps along the road to 
mapping out these fundamental quantities are already being taken at Jefferson 
Lab and HERMES. 



5 Conclusion 

Thanks to recent advances in accelerator technology that have enabled precise 
data to be collected at the world's particle accelerators, we have been able to 
probe the fascinating inner structure of the nucleon with unprecedented clarity. 
Though much has been learned from inclusive DIS experiments, future analyses 
of nucleon structure will focus more on semi-inclusive reactions, which will 
enable the spin and flavor composition of protons and neutrons to be resolved 
with greater precision. Furthermore, there is a growing appreciation of the 
need to understand the common underlying physics revealed through a range 
of observables, from elastic form factors to deep-inelastic structure functions. 

Some of the most exciting recent developments have been in the study of 
the non-perturbative structure of the proton sea through asymmetries in sea 
quark distributions, which illustrate the relevance of chiral symmetry breaking 
in QCD even at high energies. Important breakthroughs in our understand- 
ing of the proton spin have opened the way to accessing for the first time 
information about the full helicity and orbital momentum distributions in the 
nucleon. In addition, perhaps longest overdue is the need to determine the va- 
lence quark distributions in the region of large Bjorken-x, which should settle 
the long-standing puzzle of the precise x — > 1 behavior of structure functions 
and shed light on the mechanisms of spin-flavor symmetry breaking in the 
nucleon. 

To make the inroads necessary to achieve a deeper understanding of these 
issues will require full utilization of the high luminosities and machine duty 
factors at modern accelerator facilities such as Jefferson Lab. We can anticipate 
the new generations of experiments to reveal much more of the intriguing world 
of subnucleon dynamics. 
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